Langevin dynamics of the deconfinement transition for pure gauge theory 
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We investigate the effects of dissipation in the deconfinement transition for pure SU{2) and SU{3) gauge 
theories. Using an effective theory for the order parameter, we study its Langevin evolution numerically. Noise 
effects are included for the case of SU{2). We find that both dissipation and noise have dramatic effects on the 
spinodal decomposition of the order parameter and delay considerably its thermalization. For SU{3) the effects 
of dissipation are even larger than for SU{2). 
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I. INTRODUCTION 

Recent results from lattice QCD ^, corroborated by exper- 
imental data from BNL-RHIC f^], indicate that strongly in- 
teracting matter under extreme conditions of temperature and 
pressure undergoes a phase transition to a deconfined plasma. 
Such extreme conditions are believed to have happened in the 
early universe, and might also be found in the core of neutron 
stars |3]. 

The process of phase conversion during the deconfinement 
transition can occur in different ways. For a pure gauge 
SU{N) theory, the trace of the Polyakov loop provides a well- 
defined order parameter |4, 5, 6], and one can construct an 
effective Landau-Ginzburg field theory based on this quantity 
QIal- The effective potential for T <<Tc has only one min- 
imum, at zero, where the whole system is localized. With the 
increase of the temperature new minima appear (A^ minima for 
Z{N), the center of 5'L'^ {N)). At the critical temperature, Tc, all 
the minima are degenerate, and above Tc the new minima be- 
come the true vacuum states of the theory, so that the system 
starts to decay. In the case of SU{3), whithin a range of tem- 
peratures close to Tc there is a small barrier, and the process 
of phase conversion will be guided by bubble nucleation. For 
larger T, the baiTier disappears and the system explodes in the 
process of spinodal decomposition. For SU{2), the transition 
is second-order, and there is never a barrier to overcome. 

In this paper, we consider pure SU{2) and SU{3) gauge 
theories, without dynamical quarks, that are rapidly driven to 
very high temperatures, well above Tc, and decay to the decon- 
fined phase via spinodal decomposition. We are particularly 
interested in the effect of noise and dissipation on the time 
scales involved in this "decay process". In what follows, we 
adopt an effective model proposed in Ref. H for the order 
parameter and the effective potential. Numerical calculations 
for the evolution of the order parameter are performed on a 
lattice, using a local Langevin equation. 
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The paper is organized as follows. Section II briefly de- 
scribes the effective model for the order parameter. In Sec- 
tion III, we consider the Langevin evolution, discussing how 
to fix the dissipation coefficient from lattice simulations, and 
present our results for SU{2) and SU{3). Section IV contains 
our final remarks. 



II. THE EFFECTIVE MODEL 

The model proposed in | §i] intends to provide a better rep- 
resentation of lattice results for the gluon plasma equation of 
state as compared to the usual bag model. It is obtained com- 
bining a few phenomenological inputs with Z{N) symmetry 
and some known features of the perturbative equation of state. 
In particular, in a temperature range going from the deconfine- 
ment temperature Tj to 5T^i the model gives reasonable results 
and exhibits a thermodynamic behaviour that is coherent with 
data obtained from lattice simulations. 

In this approach, thermodynamic properties are determined 
by functions of the Polyakov loop, defined in Euclidean finite 



temperature gauge theories as 



qop. 
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where 'T denotes Euclidean time ordering, g is the gauge 
coupling constant and Aq is the time component of the vec- 
tor potential. We work with SU{2) and SU{3), representing 
the color degrees of freedom. Consequently, we have a Z{N) 
symmetry for the case of pure gauge theories that is sponta- 
neously broken. It would be explicitly broken in the presence 
of quarks. 

Working in the imaginary time framework, we have bosonic 
fields being periodic and fermionic fields being antiperiodic in 
the imaginary time t: 



A^(x,P) = +A^ix,0) , qix,^) = -q{x,0) 



(2) 



Any gauge transformation periodic in i respects these bound- 
ary conditions. However, as demonstrated by 't Hooft |5], 
one can consider more general gauge transformations which 



are only periodic up to the center of the group: Q.{x,^) — 

Q.C , n(x,o) = 1. 

Color adjoint fields are invariant under these transforma- 
tions, while those in the fundamental representation are not: 

A"(f,p) = Q.lA,{x,^)a, = A^(x,p) = +A,{x,0) , (3) 



^"(x,P)=nJ^(x,P)^~^(x,0) 



(4) 



Consequently, pure gauge theories have a global Z{N) sym- 
metry, which is spoiled by the addition of dynamical quarks. 

Thus the action is invariant under Z{N) transformations, 
but {TtfP{x)) is not. Symmetry requires z {Ti fP{x)) = 
{TtfP{x)), which implies (Tr/? P(x)) = 0. When the phase 
transition occurs, the Z{N) symmetry is spontaneously bro- 
ken and (Trf f (x)) assumes a non-vanishing value. So, one 
can use it as an order parameter for the transition and it is 
possible to write an effective theory for its dynamics. 

The efective theory of Ref. |^ is based on a mean field 
treatment in which the Polyakov loops are constant through- 
out the space and the free energy is a function of its eigenval- 
ues. A perturbative calculation of the free energy of gluons as 
a function of the Polyakov loop eigenvalues yields 
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where 9 is defined through the eigenvalues of the Polyakov 
loop, Pjif = exp(/9;) hjk, and AG^^^ = Qj — 9^, which reduces 
to the usual blackbody formula in the case Aq = 0. 

This expression for the free energy of gluons propagat- 
ing in the background of Polyakov loops predicts a gas of 
gluons that is always in the deconfined phase, with no in- 
dication that higher-order corrections will modify this re- 
sult. This can be modified by the introduction of a mass 
scale into /. This mass scale will determine the deconfine- 
ment temperature 7j and is introduced in a phenomenological 
way. One ends up with the same expression of the pertur- 
bative calculation, but now w^. = \/W+M^. Parametrizing 
the Polyakov loop and representing the diagonal matrix as 
rfia^[exp(/(|)^/2,...,/(|)i,— /(|)i, ...,— /(|)Ar/2], it is possible to ex- 
tract an effective potential as a function of (|). For SU{2) one 
has only (|)i = (|) and (|)_i — — (|), and the effective potential can 
be written as 



V = -- 



15 
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.... ... -TT). (6) 

4 Tf^ 

Notice that there is a symmetry (|) <-^ 7i — (|) associated with 
the Z(2) invariance. It is convenient to write this equation in 
terms of a new variable x|/ = 1 — isfn/l to make this symmetry 
more evident. One obtains: 



V = -- 
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where \|/ = represents confinement. It is interesting to con- 
nect the \\t used here and the trace of the Polyakov loop, 
used in Ref. |7] as the order parameter For the diagona- 
lyzed matrix we have the trace, according to our parametriza- 
tion, as e"^ + e^'^, namely TrL = 2cos((|)), or, as defined 
above, TrL = 2cos(7t(l — V(/)/2). So, when x|/ = we have 
TrL = 0, which represents confinement, and when \\t^l, then 
TrL -^ 1, representing the deconfined state. 

The phase transition in this case is second order, as 
expected |9]. The value of M can be determined from 
the deconfinement temperature through the relation L^ — 
{3/2yi^M/% « 0.38985M, so that it is possible to extract the 
deconfining temperature from the lattice and then fix M. The 
minimum of the potential occurs for 
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For5L'^(3) there are three eigenvalues: (|)i = (|), and (|)_i 
-isf. The potential assumes the form: 

o— 2 T^S 

V = -r3-- + _[8^2(^_^)2^^2(^_27r)2] 
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Again, it is useful to rewrite the potential in terms of a new 
variable \|/ = 2jt/3 — (|), so that one obtains 

405 \2%^ ^ 



37t ^ 27t2 ^ 

Now, M and Tj are related as follows: 
9 



Td = ^VIOM K. 0.45296M , 
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In this case, TrL = e'^ + \+e '^ and the connection with \|/ 



becomes TrL = | cos 
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Immediately above the critical temperature, the SU{y) po- 
tential presents a barrier between the old and the new vacua. 
This barrier, however, is very small and quickly disappears 
with the increasing of the temperature. One should notice that 
above 2Tc the changes in the potential are negligible. 



III. LANGEVIN EVOLUTION 

Let us now consider the real-time evolution of the order 
parameter for the breakdown of Z{N). We assume the system 
to be characterized by a coarse-grained free energy 



»,n^/ 



= / d'x 



^(v^)2 + y,^^(c^,r) 



(13) 



T = 216McV 
T = 270 MeV / 
T = 378 MeV / 
T = 500MeV / 




FIG. 1: Effective potential for SU(3) for different values of the tem- 
perature. 



where Veff{^, T) is the effective potential obtained in the last 
section, and B = K^T/g^ for SU{2) and B = AT/g^ for SU{3). 
The time evolution of the order parameter and its approach to 
equilibrium will be dictated by the following Langevin equa- 
tion 
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where g is the QCD couphng constant, and Y is the dissipation 
coefficient, which is usually taken to be a function of temper- 
ature only, r = r(r). The function ^ is a stochastic noise 
assumed to be gaussian and white so that (^(x,f)) = and 

(^(x,r)^(x',f')) = 2r5(x-jc')5(f -f')- The noise and dissi- 
pation terms are originated from thermal and quantum fluctu- 
ations resulting either from self-interactions of the Polyakov 
loop field or from the coupling to diferent fields (such as chi- 
ral fields). The case with only first-order time derivative was 
considered in Ref . llOll . 

This description is admittedly very simplified. A more 
complete analysis should consider different contributions of 
noise and dissipation terms and memory kernels instead of 
simple Markovian terms proportional to the first time derivate 
of the field |11, 12]. In general, one obtains a complicated 
dissipation kernel that simplifies to a multiplicative dissipa- 
tion term which depends quadratically on the amplitude of the 
field as Fi [T) Vf/^ {x,t)L{x,t) where Fi is determined by imag- 
inary terms of the effective action for \(/ and depends weakly 
(logaritmically) on the couplings. The fluctuation-dissipation 
theorem implies, then, that the noise term will also contain 
a multiplicative contribution of the form \|/(i',f)^(x,f ), and be 
in general non-Markovian. The white noise limit is reobtained 
only for very high temperatures. 

For the SU{2) case we have fixed F in the following way. 
We have used pure-gauge Euclidean lattice Monte Carlo sim- 
ulations in the line discussed in Ref. j^. In this approach, 
spinodal decomposition is obtained on the lattice performing 
local heat-bath updates of gauge field configurations at p = 
4/^2 = 3, after thermalizing the lattice at p = 4/^^ = 2. The 
critical value of p for deconfinement is found to be p^ ^ 2.3. 



F is then extracted by comparing the short-time exponential 
growth of the correlation function {L{k,t)L{~k,t)) predicted 
by the lattice simulations |13|| and the Langevin description, 
assuming of course that both dynamics are the same. Making 
this comparison for the lowest lattice momentum mode, it is 
found that F = 7.6 x 10^ T^ //j, where /j is a time scale relat- 
ing Monte Carlo time and real time. Assuming that typical 
thermalization times are of the order of a few fm/c, we obtain 
F- lO^fm-2. 
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FIG. 2: Volume average of the SU{2) order parameter normalized to 
the \|/o > minimum of the bare effective potential. 

In our numerical calculations we solve Eq. J14t on a cubic 
spacelike lattice with 64^ sites under periodic boundary con- 
ditions. We use a semi-implicit finite-diference scheme for the 
time evolution and a finite-difference Fast Fourier Transform 
for the spatial dependence ] 14]. For SU{2), the critical tem- 
perature is Td = 302 MeV ] 15] and we obtain M = 775 MeV. 
We took the average of several realizations with random ini- 
tial configurations around i|/ ^^ 0. We consider the time depen- 
dence of the volume average of \\t 



(v) = ^Ev'>*(0 



N^ 



(15) 
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where A^ is the number of lattice sites in each spatial direction, 
and i,j,k = !,■■■ ,N are the lattice sites. In Fig. |2lwe plot 
(i|/)/i|/o, where \|/o > is the positive minimum of the bare 
effective potential, for three situations: no dissipation and no 
noise (dotted curve), no noise (dashed curve) and full solu- 
tion (solid curve). When considering noise, we have added 
the appropriate counterterms to make the equilibrium solu- 
tion independent of the lattice spacing 1 16] . All curves are for 
r = 6.6 71/. 

Clearly seen in Fig. |2l is the large effect of dissipation, 
which delays the rapid exponential growth of the order param- 
eter due to spinodal decomposition. The retardation seen here 
for the deconfinement transition is substantially larger than the 
corresponding delay seen for the chiral condensate evolution 
in Ref. ] 17]. The effect of noise is also in the direction of de- 
laying equilibration, as expected. Also expected, and clearly 



shown in Fig.|2] is the effect of noise in the equilibrium value 
of i|/ which is larger than i|/o. 
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FIG. 3: Volume average of the 5(/(3) order parameter normalized to 
the \(fo > minimum of the bare effective potential. 

For 5C/(3) one has Tj = 263 MeV (H, so that M ^ 
580 MeV. The results of our simulations at T = 6.6 Tj are 
shown in Fig.|3] Here we are using the same lattice and same 
dissipation F as for SU{2). As seen in this figure, the effect 
of dissipation is even more dramatic than for SU{2), with the 
proviso of course that we are using the value of F extracted 
from SU{2) lattice simulations. As mentioned earlier, im- 
mediately above the critical temperature the SU{3) potential 
presents a barrier between a local minimum and an absolute 
minimum. However, this barrier has no effect on the delay 
seen in Fig.|3] since our simulations are done for high temper- 
atures, T ^ 2Td. We have not investigated the effect of noise 
in this case because the appropriate renormalization counter- 
terms for an effective potential with a first-order transition are 
not yet available 1161 . 



rV. SUMMARY AND OUTLOOK 

We have investigated the effects of dissipation and noise 
in the deconfinement transition of SU{2) and SU{'i) pure 
gauge theories. We have used the effective model proposed 
in Ref. 18], which combines phenomenological inputs with 
Z{N) symmetry and some known features of the perturbative 
equation of state. The model provides a reasonable represen- 
tation of lattice results for the pure-gluon plasma equation of 
state in the temperature range between T^ and 57"^ . We have 
performed numerical simulations for the evolution of the or- 
der parameter on a spatial cubic lattice using a local Langevin 
equation. We find that both dissipation and noise have dra- 
matic effects on the spinodal decomposition of the SU{2) or- 
der parameter, delaying considerably its thermalization. Dis- 
sipation effects are even larger for 5t/(3). 

The present work must be improved in several aspects. Per- 
haps the most important one is in the method used to extract 
the dissipation coefficient F PIS']. This was done using Eu- 
clidean lattice Monte Carlo simulations, in which spinodal de- 
composition of the order parameter is obtained performing lo- 
cal heat-bath updates of gauge field configurations above the 
deconfinement temperature. One of the major uncertainties 
in this approach is the relation between Monte Carlo updates 
and real time. Another source of uncertainties comes from a 
richer structure of noise and dissipation terms, including an 
evaluation of memory kernels. It is widely known that, in 
general, quantum corrections lead to complicated dissipation 
kernels that only in very special situations simplify to an addi- 
tive noise term as used here. These issues will be considered 
in a future publication 1 19] . 
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